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Abstract 

In the hadrodynamics (QHD) there are two methods to take ac- 
count of the contribution of negative-energy states in the relativistic 
random phase approximation (RRPA). Dawson and Furnstahl made 
the ansatz that the Dirac sea were empty, while according to the Dirac 
hole theory the sea should be fully occupied. The two methods seem 
contradictory. Their close relationship and compatibility are explored 
and in particular the question of the ground-state (GS) instability 
resulting from Dawson-Furnstanhl's ansatz is discussed. 
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1 Introduction 



In the framework of the conventional quantum hadrodynamics (CQHD), the 
relativistic mean-field theory (RMFT) based on effective Lagrangians and the 
no-sea approximation has achieved a remarkable success for the calculation 
of nuclear ground states (referred to as MFGS i.e. mean-field GS) [1]. Since 
the parameters are adjusted to fit the experimental data, one expects that 
a large part of vacuum contributions and some other neglected effects are 
already taken into account phenomenologically [1]. A formal justification 
of the above conjecture has been given by Furnstahl, Pickarewicz and Serot 
[2] on the basis of the effective field theory (EFT), where a more elaborate 
effective Lagrangian with non-renormalizable interactions (modern QHD [3]) 
has been considered. For the calculation of excited states in the relativistic 
random phase approximation (RRPA), Dawson and Furnstahl (D-F) pointed 
out long ago [4] that in order to separate out the spurious J 71 " = 1~ state 
and to preserve current conservation, it is of vital importance to take the 
negative-energy (NE) states into account, because only then will the Dirac 
single-particle basis become complete. They suggested that based on the 
above MFGS one may assume the NE sea is empty and thus besides the 
positive-energy (PE) particle-hole pairs one should further consider pairs 
formed from a particle at one of the NE states and a hole in the occupied PE 
states (referred to as ah pairs as in Ref. [5]). Ma et al. [6] found that the D-F 
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method can improve the calculation of nuclear giant resonances significantly. 
Recently, Ring at al. [5] have made a detailed study of the effects of the ah 
pairs and found that they are indeed important. However, as is wellknown [7, 
8], an empty NE sea cannot avoid the problem of spontaneous radiation, i.e. 
of keeping the valence nucleons from jumping into NE states through emitting 
mesons and photons, if the corresponding interaction terms are included in 
the Lagrangian. As numerous NE states are available, the MFGS, even 
though successful in all the other aspects, will not be stable and will disappear 
almost instantly. For this reason we would like to emphasize that it is still 
important to follow Dirac's original proposal of a filled NE sea. Evidently, a 
filled NE sea can equally fulfill all the above requirements pointed out by D-F. 
Besides, if any effect due to the filled sea has been encoded into the effective 
Lagrangian by fitting the parameters to experimental data, in the calculation 
it should and can be neglected in order to avoid redundancy (see section 2). 
Clearly, for such a purpose there is no need to drive the sea particles out of the 
sea. Since the calculation on the basis of a filled sea will be as simple as one 
with an empty sea if the vacuum loops are neglected, it shows that Dirac's 
postulate is a preferable choice for both CQHD and MQHD, because it cannot 
only prohibit the valence nucleons from tumbling into the NE sea, but also 
take account of the existence of antiparticles in a natural way. In addition, it 
also provides a greater possibility to treat the vacuum contribution effectively, 
namely if it is already included in the parameters of the effective Lagrangian 
(EL), we may simply neglect it in the calculation; however, if a part of it 
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has been left out owing to a somewhat excessive truncation, we may then 
either add additional terms and parameters to EL or calculate the part as a 
contribution from the sea. It is known that the results of the one-N (nucleon) 
-loop and virtual NN pair calculations based on structureless nucleons are 
questionable and earlier calculations with CQHD met difficulties of large 
effects from loop integrals caused by the vacuum contribution. It seems 
[9] that the above drawbacks can be overcome by a careful study of vertex 
corrections and relevant higher-order effects. 

The relation between the effective D-F method and Dirac's hole theory 
will be studied in some detail in section 2. In section 3 numerical results 
calculated in RRPA are presented and discussed, where phenomenological 
vertex corrections are also considered to make a preliminary study of the 
effect due to the nucleon compositeness. Finally a summary is given in the 
last section. 

2 Relation between D-F's effective method 
and Dirac's hole theory 

In order to expose the relation between D-F's effective method and Dirac's 
hole theory more clearly and in simpler terms, we shall restrict our discussion 
to symmetric nuclear matter and Walecka's er — uj model. The Lagrangian is 
written as 



L — — 
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where W^ v = d^Uy — d v uo^ and ip, a, on denote the nucleon, a-meson and lo- 
meson fields, respectively. As is wellknown [1] , through fitting the parameters 
to nuclear properties, the a — uj model, though simple, has achieved consid- 
erable success in the description of many bulk and single-particle properties 
of nuclei. Thus we may assume that a part of vacuum contributions and 
some higher-order effects have been included in the Lagrangian L, though 
the inclusion cannot be so complete as what the more elaborate effective 
Lagrangian may provide. We would like to show that if the vacuum contri- 
bution from Dirac's filled NE sea is neglected, the Dirac results reduce almost 
to what one obtains according to the D-F ansatz. In order to explore the 
relationship between the two methods we may apply them to calculate the 
same set of physical quantities and compare their results. For this purpose 
we shall consider the correlation function 

C(A, B; x u x 2 ) = (T[A(x 1 )B(x 2 )}} - (A(x 1 )) (B(x 2 )} , (2) 

as a typical example. In Eq. (2) (O) = (^o |0|^o) > K(x) = ip(x)T K ip(x) 
(K = A or B), ip(x) is the nucleon field operator, whereas operator T K is field- 
and time- independent. For instance, for the isovector multipole operator we 
have T K = •y 4 Tsr x Yx fJ ,(0,(p) and for the bilinear Dirac current and scalar 
density F K — 7 (k — /i — 1, 2, 3, 4) and F K — 1 (k — 5), respectively, 
etc. For simplicity we shall further assume that T K contains no differential 
operators. In the lowest order approximation Eq. (2) has the form 

C°(A,B;x 1 ,x 2 ) = -Tr[T A (^ 1 )G°(x 1 - x 2 )T B (x 2 )G°(x 2 - x,)], (3) 
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where G°(x) is the relativistic Hartree approximation to the nucleoli propa- 
gator 

G aP {x = xi- x 2 ) = (T[-0 Q (xi)^(^ 2 )]) , (4) 

and x = x li = (x, ix ) with x = t. The Fourier transform of Eq. (3) is given 
by 



C°(A,B;kk') 



JJ d 4 x 1 d 4 X2e 
—2n5(ko — k 



ikxi+ik ' X2 C°(A,B;x 1 ,x 2 ) 
d A q d 3 p 



Tr[r A (k-k'-p) 



(2tt) 4 (2tt) 3 

xG°(k' + q + p, k' + go)rB(p)G°(q, q )}. (5a) 

If T K is further independent of x (indicated by k taking a small letter), we 
have r a (p) = (27r) 3 £ 3 (p)r a and Eq. (5a) reduces to 

d 4 q 



C°(a,b;kk>) = (2rr)Vf ' ! /,- //) ; ^ 
= (2tt) 4 5 (4) (A;- k')C°(a,b;k) 



Tr[r a G°(q + k)r b G°(q)] 



(5b) 



One observes that except for the 5-function and a constant factor, Eq. (5b) 
is just the expression for the polarization tensor in the a — u model, if F K 
(k — 1 to 5) = {7^, 1}. Since we are considering a nuclear matter whose GS 
is given by RMFT, G°(p) = G° F (p) + G° D (p) can be written as 



g° f (p) = -(w + i*r) 



P 



1 



p 2 + M *2 

_ { 0(E P 



+ m5{j> 2 + M* 2 ) 

e(-E T 



2\E p \ p -E p + ie p -E p -ie 



(6a) 



} (6b) 



ITT 



gUp) = trip* + iM*)—e{ PQ )9{k F - |p|)5(p - e p ), 



(6c) 
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where E p = ±[p 2 + M* 2 ] 1/2 , M* is the effective mass, hp denotes the Fermi 
momentum, F the Feynman part, D the density-dependent part and P the 
principle value, we may rewrite Eq. (5) as 

C°(A, B; kk') = C° FF (A, B; kk') + C° m (A, B; kk') + C° DD (A, B; kk'), (7a) 

C^(A, B; kk') = C° DF (A, B; kk') + C° FD (A, B; kk'). (7b) 

On the right-hand side of Eq. (7) the subscript D or F indicates G° in Eq. 
(5) is G° D or G° F and the three parts will be referred to as FF-, m- and De- 
part, respectively. It is wellknown that ip(x) = e lHt ip(x.,0)e~ lHt and ?/>(x, 0) 
can be expanded in the form 

4 j3 

V(x, o) = E / j^<P r Mp r y p ^ (8) 

where we note c{pr) = d + (—p,r — 2), u(pr) = v(—p,r — 2) if r = 3, 4 and 
v(pr) (r = 1, 2) denotes a NE spinor with momentum —p. Substituting Eq. 
(8) into Eq. (4), we easily find [10] that the Fourier transform of G af 3(x) can 
be represented in the form 

4 

Gapip) = X) u a (pr)iG(pr,ps;p )uf3(ps), (9a) 

r,s=l 

G(pr,ps;p ) = - dt e tpot {T[c{prM)c + {psM)]) ■ (9b) 

Hereafter, we shall use a tilde to indicate the result obtained according to 
the D-F method. D-F's ansatz [4] implies 



G°(pr,ps;p ) = 6 



9(E P - E F ) 6(-E p ) , 9(E F - E p )9(E r 



p - E p + it p - E p + it p - E p - it 
S rs G°(p). (10a) 



Inserting Eq. (10a) in Eq. (9a), we obtain G°(p) = G° F (p) + G° D {p) , G° D (p) = 
G° D {p) and 

^,. (7A+un ^{_i^___S^_}. (10b) 

It is seen that the only difference between G°(p) and G°(p) consists in their 
F part. In G° F (p) , just as G°(p) in Eq. (10a), its NE pole is now shifted to 
the lower-half plane and there are no other changes from G F (p) to G° F (p). 
Defining AG°(p) = G°(p) - G°(p), we have From Eqs. (6) and (10) 

ITT 

AG°(p) = AG F (p) = (7^ + tM*)—5( Po - E p )9(-E p ). (11) 

Since all the poles of G° F (p) are in the lower-half plane, substituting Eq. 
(10) in Eq. (5), as pointed out in [4], we get C FF (A,B; kk') = which also 
follows from the fact that here there are no antiparticles. Besides, we find 
C° DD (A, B; kk') = C° DD (A, B; kk') and 

C° m (A, B; kk') = C° m (A, B; kk') + AC° m , (12a) 

AC = C° DAF (A, B- kk') + C° AFD (A, B; kk'). (12b) 

In Eq. (12b) the subscript AF means that we substitute AG° F (p) for G° F (p) 
in Eq. (7b). Comparing Eq. (12) with Eq. (7), we have 

AC = C° - C° = AC° m - C° FF , (12c) 

i.e. the ansatz of empty NE sea is equivalent to the assertion that C FF in C° 
is replaced by AC° t . Thus, if the vacuum contribution C FF can be neglected, 
the difference between the two results consists solely of AC^. Note Eqs. (7) 



and (12) show why the D-F method can give a better result than the no-sea 
approximation, because the latter only takes account of C^ D and a part of 
given by the first term in Eq. (6b) or (10b), while the former has further 
correctly considered the full with even a correction AC^ (see below). In 
order to expose the physical implication of Eq. (12a) more clearly, consider, 
for instance, Eq. (5b) with T a = 7^ and Tj, = j x . Substituting Eq. (6a) in 
Eq. (7b), and using a prefix 5 to denote the contribution of the second term 
of Eq. (6a), we obtain 

5C° m (r), A; k) = ReC ^, A; k) = M 1 (r ] , A; k) + M 2 ( V , A; k), (13) 

Mi(rj,\;k) = / A^^{^^^o + * q - £ q+k ) 

+ t ^~ k) 5{k - £ q + E q _ k )}, (14a) 

+ t ^y k) S(k a -E^-E^)}, (14b) 

where Re means the real part and 

t vX (q ± k) = [2q v q x ± q v k x ± q x k v =F 6 vX (q • k)} qo=E ^. (15) 

From Eqs. (11-13), it is easily seen that 

AC° m (vA;k) = -2M 2 (r ] ,X;k), (16) 

Re&Jj], A; k) = ReC° m ( V , A; k) + AC^fa, A; fc) = M 1 {r } , A; fc) - M 2 ( V , A; fc), 

(17a) 



ImC° m (r), A; k) = ImC° m ( V , A; k). (17b) 

According to Eqs. (16) and (17) one concludes that and differ from 
each other only in the sign before M 2 , because their imaginary parts are the 
same. If we require k > 0, it is seen from Eq. (14b) that the first term 
in M 2 is zero and only the second term will contribute. The (^-function in 
Eq. (14b) shows M 2 becomes effective only if k > 2M*, which is about 
1.37GeV for M* ~ 0.73M, i.e. if the excitation energy is not too high, 
we have (5° t ~ C Q m . Since M 2 = if < k < 2M*, in the energy region 
where for instance, giant multipole resonances are studied we even expect 
= C^. Clearly the above derivation also applies to the other components 
of C°(a, b; k) (see Eq. (5b)) as well as to Eq. (5a) and similar results obtain. 
Sofar we have only considered the lowest order approximation to Eq. (2). 
In order to gain some idea about the accumulative effect, we have further 
made a RRPA calculation for C(a, b; x 1: x 2 ) with C° given by Eq. (5b). The 
Dyson-Schwinger equation for the full scalar-vector meson propagator in the 
a — oo model can be written in the form [1 1] 

Aab(k) = A° ab (k) + Al(k)U° cd (k)A db (k), (18) 

where the polarization tensor n° 6 (fc) is related to C°(a, b; k) in Eq. (5b) with 
T K = {7 , 1} for k = a or b as follows: 




K U U\ = ( 9 2 v C°(fi,iy;k) -ig v9s C\^h-k) \ 

Ul n° 5 ) { - igs g v C°(5,u;k) -g 2 s C°(5,5; k) )' 



(19) 
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If we substitute T K = {ig v lfj,, 9s} for T K and use an overhead bar to indicate 

an expression obtained in this way, we have C (a,b;k) = -Tl° ab (k). In Eq. 

/ a°5 \ 

(18) the free meson propagator is given by = I ^ '" ^ v ^ J, where 
A°(/c) = —i[k 2 +mf.—ie}~ 1 for k = v or s, and m K is determined by [k 2 +m 2 + 
Hl° K (k 2 ; F)| M . =M ] fc2= _^2 = with m K denoting the physical mass, and for 
our present purpose no medium dependence of meson masses is considered. 
It will be studied in more detail elsewhere. If only the summation of ring 
diagrams is considered, as illustrated in Fig. 1, RRPA to C(A, B;kk') (see 
Eqs. (2) and (5a)) may be written as 

C(A, B; k, k') = C°(A, B; k, k')- J -^C°(A, c; k, q)A cd (q)C°(B, d; —k', -q), 

(20a) 

C°(EJ;k,q) = 27c5(k -q ) J ^Tr[r £ (k - q)G(p + q)T f G(p)}, (20b) 

where E = A or B. Eq. (18) shows the ring summation is achieved by A c ^ 
in Eq. (20). If T E is independent of x, from Eq. (18) one finds easily by 
iteration that Eq. (20) can be simplified to CjiPA(a,b; k) = —H ab (k) with 

IWA;) = K b (k) + U ac (k)A° cd (k)U db (k). (21) 

3 Numerical results 

We have solved Eq. (21) with n° 6 given by Eq. (19), where the scalar and 
vector mesons as well as their mixing are considered. The parameters are 
g 2 = 0.5263, g\ = 0.6842 (g 2 = ^ 2 /16tt 2 ), whereas M* = 0.73M, M = 939, 
m s = 550, m v = 783 (all in MeV), and k F = 1.3/m -1 [1]. 
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To illustrate the effects of C FF and A(7^ (see Eqs. (7) and (12)) as well 
as the difference in results between the Dirac and D-F theories, we shall, 
as an example, discuss (7(4) = ite[C(4, 4; k)] = g~ 2 Re[Tlu(k)], which is 
nonnegative and closely related to the longitudinal response function. Fig. 2 
shows the dependence of (7°(4) and (7°(4) on the energy transfer k . Indeed, 
in the small ko region, we have (7°(4) = C' (4), because both C FF (A) and 
A (7^ (4) will not be zero only if k > 2M*. Thus, to the lowest order, the 
two methods yield the same results in case k < 2M*. In the large k region 
since C° DD {4) = C° DD {4) = (see Eq. (A. 26) in Ref. [13]) and according 
to Eq. (14a) Mi (4, 4; k) = for time-like k^ [11], one gets <7°(4) = (7^(4), 
C^(4) = -(7^(4) and <7°(4) = (7^(4) + C° FF (A). From Fig. 2 it is seen that 
the sign of (7^(4) is correct, while (7^(4) has a wrong sign. However, in the 
Dirac theory the relevant effect is represented by (7° (4) whose sign is again 
correct owing to the contribution of C FF (A). It is interesting to find that 
A (7^ included in C Q rn in Eq. (12a) actually means a correction rather than a 
defect, though its effect is not very important. For the RRPA calculation we 
shall consider three cases: (1) Both T a and T b in the Eq. (5) are bare vertices. 

(2) One of F K (k — a or b) is a dressed vertex T K , which is phenomenologically 
taken as 

f K = T Kj ^-^T K F(k^A). (22) 

(3) Only one T K in C FF is replaced by T K . For simplicity, the case of dressed 
vertex at each end of the loop will not be considered. The detailed formula 
for II° 6 and the relevant renormalization procedure are known for all the 
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above three cases [11-14], thus they will not be written down here. We shall 
use a subscript r\ = 1 to 3 to indicate the case which the calculation refers 
to. In Fig. 3 the calculated RRPA results for cases (1) and (2) are shown, 
where we have chosen A = l.lGeV for case (2). Since according to Eqs. (18, 
19) and (21) the real as well as the imaginary part of n° d (/c) contributes, it 
is seen that there is a significant difference between Ci(4) and Ci(4) even in 
the small k region. Hereafter we shall restrict our discussion to this region, 
as the parameters in the effective Lagrangian are fitted to nuclear ground- 
state properties. Fig. 3 shows that the difference between Ci(4) and C 2 (4) is 
noticeably larger than that between Ci(4) and 6*2(4). It shows that the form 
factor given in Eq. (22) is rather ineffective in the Dirac case. However, if we 
assume that the effect of nucleon compositeness has already been encoded 
in the effective Lagrangian, we may simply set F(k 2 ,A) = 1 i.e. then there 
would be no need to consider C 2 (4) and (5 2 (4) in order to avoid redundancy. 
From Eqs. (7) and (19), it is easily seen that the large difference between 
Ci(4) and Ci(4) is caused by C FF calculated with bare vertices. For this 
reason we have further considered case (3), which is represented graphically 
in Fig. 4, where the graphs of case (1) are drawn for comparison. We shall 
use C 3 (4,o;) to denote the results calculated with T K (a) = aT K F(k 2 , A e ), 
where a (— to 1) is a percentage factor and A e = l.lGeV. If a — 1, we 
have C 3 (4, 1) ~ C 2 (4) ~ C*i(4), i.e. the result of Dirac's case. On the other 
hand, if a — 0, we have T k (a) = 0, thus C 3 (4, 0) = C 1 (A) as shown in Fig. 
4, i.e. if C FF or the vacuum contribution is neglected, the Dirac and D-F 
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results are the same, because the effect of AC^ is insignificant. Note Fig. 
4(b) shows that in the case of k a = 330MeV C 3 (4, a) varies with a unusually. 
It first rises and then falls down finally to C 2 (4) ~ Ci(4) at a — 1. Since 
Ci(4) ~ C 2 (4), this suggests that (7 2 (4, a) calculated with T K replaced by 
T K (a) may display such a similar behavior. This is indeed so as shown in 
Fig. 3(b) where (7 2 (4, 1) = 6*2(4). Clearly, the postulate of filled sea also 
gives a greater possibility to realize the experimental data than the ansatz 
of empty sea. In fact, if for instance, C 3 (4,a; = 0.1) fits more closely to 
the experimental data, it means that the effective Lagrangian has not yet 
included the vacuum effect and nucleon compositeness completely and there 
is no need to add additional terms and parameters to the Lagrangian, as they 
can be taken into account by means of Dirac's method in the simple way. In 
the above, we have argued and demonstrated if the effective Lagrangian has 
encoded the vacuum contribution and some other higher-order effects (for 
instance, nucleon compositeness) in its parameters, we may simply neglect 
these contributions in the calculation to avoid redundancy without assuming 
that the NE sea is empty and nucleon is structureless. 

4 Summary and conclusions 

To compare the Dirac and D-F methods, we have considered the calculation 
of correlation function as a typical example. As shown in section 2, if the 
vacuum contribution can be neglected, the Dirac and D-F results are the 
same in the small ko region not only to the lowest order but also in RRPA, 
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because according to Eq. (16) the effect of AC^ is insignificant in the model 
studied. Thus, in such cases, the D-F ansatz will yield a good effective 
method of calculation, though it is unnecessary, because the Dirac theory 
gives the same result and the calculation will be no more complicated if 
the vacuum effect and some higher order effects can be neglected. Since 
the Dirac theory can identify all the higher order effects unambiguously, to 
neglect them we only need to neglect the relevant terms, as displayed in 
section 2. Actually, no-sea approximation does not imply that the sea must 
be empty. It only means that in the calculation the vacuum contribution 
can be neglected, because it has already been encoded in the parameters of 
the Lagrangian. As is well known, there is no way to avoid the spontaneous 
radiation which will cause MFGS to disappear instantly if the sea is empty. 
We have illustrated in sections (2) and (3) why we have emphasized in the 
introduction that it is not only important, but also advantageous to follow 
Dirac's original postulate of filled sea. Clearly if the effect of AC^ is not 
small, the D-F and Dirac results may differ, even if the vacuum contribution 
can be neglected, thus the difference may then serve as a discrimination 
between the two methods. This is a problem under further study. 

This work is supported in part by the National Natural Science Founda- 
tion of China and the Foundation of Chinese Education Ministry. 
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Figure Captions 

Fig. 1: Graphical representation of RRPA for the calculation of the 
correlation function C(A, B; Xi,x 2 )- 

Fig. 2: Real part of the lowest order approximation to the correlation 
function C(4,4; k): C°(4) = Re[C°(4,4; k)}. In the small k region the three 
curves coincide. 

Fig. 3: Curves represent C(4) = Re[C(A,A; k)] in RRPA. The solid 
curve denotes Ci(4), the thin one Ci(4), whereas the heavily dashed one is 
C 2 (4), and the dashed one (5 2 (4, a) where <5 2 (4, 1) = C 2 (4). (both with A = 
l.lGeV). (a) momentum transfer k a = |k| = 560MeV, (b) A; a = 330Me\/. 

Fig. 4: Effects of vacuum contribution and nucleon compositeness in 
RRPA. The solid and thin curve again reprint Ci(4) and Ci(4), respectively. 
The dashed curves attached with the value of a (a — 0.3, 0.6 and 1) denote 
the respective C 3 (4). (a) A; a = 560Me\/, (b) k a = 330MeV. 
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